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ABSTRACT 

This  report  describes  a  "lumped  mass"  computer  simulation  of 
a  trailing  wire  system  towed  by  an  orbiting  aircraft.  The  USAF  U-10 
aircraft  flight  conditions  were  used  as  applied  to  three  flightpaths: 
(1)  Constant  altitude  circular,  (2)  Constant  altitude  elliptical, 
(3)  Non-constant  altitude  circular  orbit.  These  flightpaths  were 
flown  in  no  wind  and  30  knot  wind  conditions. 

Wire  lengths  of  4000  ft  and  8000  ft  were  used  where  the  free 
enti  was  attached  to  a  circular  sphere  of  25  and  100  lbs. 

The  results  from  12  orbiting  trailing  wire  configurations  are 


tabulated  and/or  plotted.  One  of  the  main  objectives  in  this  analysis 
was  to  determine  how  yo-yo  of  the  free  end  of  the  wire  in  wind  might 
be  alleviated  or  at  least  minimized. 
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INTRODUCTION 

An  orbiting  trailing  wire  system  consists  of  an  aircraft  flying 
in  a  closed  path  and  towing  a  long  cable.  The  shape  and  size  of  the 
orbit  of  the  free  end  of  the  cable  are  of  prime  importance.  More 
specifically,  the  yo-yo,  or  vertical  oscillation  of  the  free  end, 
must  be  controlled  within  certain  limits  if  the  trailing  wire  system 
is  to  be  useful. 

This  paper  describes  a  computer  simulation  of  the  system,  in 
which  the  cable  is  modeled  by  a  series  of  lumped  masses  connected 
together  by  linear,  massless  springs.  The  coupled,  second  order  dif¬ 
ferential  equations  of  motion  of  the  lumped  masses  are.  programmed  onto 
a  digital  computer.  This  paper  uses  the  computer  simulation  to  study 
the  steady  state  motion  of  the  cable  with  primary  emphasis  given  to 
the  motion  of  the  endpoint  of  the  cable.  Since  there  are  only  a  finite 
number  of  mass  points  in  the  model,  errors  are  introduced  into  the 
analysis.  However,  these  errors  are  insignificant  even  if  relatively 
few  mass  points  are  taken. 

Any  orbiting  trailing  wire  system  can  be  analyzed  with  the  com¬ 
puter  simulation.  This  paper  contains  the  analysis  of  a  specific 
trailing  wire  system  referred  to  as  the  SNOWBIRD  system.  The  SNOWBIRD 
system  has  been  flight  tested  by  Air  Force  Academy  personnel  using  a 
U-10  as  the  towing  aircraft  and  a  cone  attached  to  the  end  of  the 

/ 


cable.  The  cable  in  the  computer  simulation  is  the  same  cable  that 
has  been  flight  tested.  During  the  simulation  the  aircraft  is  flown 
at  a  velocity  and  a  bank  angle  used  by  the  U-10  in  flight  tests. 

PROGRAMMING  THE  EQUATIONS  OF  MOTION 

2.1  Description  of  Lumped  Mass  Model 

The  orbiting  trailing  wire  is  modeled  by  a  series  of  lumped 
masses  connected  together  by  linear,  massless  springs.  Figure  1  is 
the  model  of  a  trailing  wire  cable  approximated  by  10  masses  and  10 
identical  linear,  massless  springs.  The  leading  end  of  the  cable,  or 
Station  1,  represents  the  towing  aircraft.  One-tenth  of  the  total 
mass  of  the  cable  is  located  at  each  of  the  10  lumped  masses.  Stations 
2  thru  11.  The  mass  of  the  drogue  is  included  in  the  final  lumped 
mass,  Station  11. 

The  equations  of  motion  for  the  trailing  wire’s  masses  are  devel¬ 
oped  in  Appendix  I  of  this  paper.  The  three  equations  of  mot  '..'I  Eor 
mass  point  i  are: 


S  S 

°i  °i 

Vi  ■  ki<1-^)(xi+rxi)-ki-i(1-  <Wi> 
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where:  «  mass  of  point 


-  spring  constant  for  segment  i 


e^  *»  stretched  spring  length 
S  «  unstretched  spring  length 


«  aerodynamic  lift 


«  aerodynamic  drag 


g  *»  acceleration  of  gravity 

V  V  \ 

in  the  lift  direction 


n  ,  i\r  ,  n„  are  direction  cosines  of  n,  the  unit  vector 
x.  1 


,  5  »  and  £  are  direction  cosines  of  5,  the  unit  vector 

zi 


V  >i 


in  the  drag  direction. 

During  each  step  of  the  integration  of  the  equations  of  motion,  the 
aircraft  is  flown  an  incremental  distance  and  the  equations  are  inte¬ 
grated  using  a  Runge-Kutta  routine.  The  displacement  of  the  aircraft 
enters  into  the  three  equations  of  motion  for  the  mass  point  immedi¬ 
ately  below  the  aircraft,  or  at  Station  2. 
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2.2  Aircraft  Orbit 


The  cowing  aircraft  is  programmed  to  fly  m  either  a  circular 
or  an  elliptical  orbit.  The  velocity  of  the  aircraft  is  constant 
with  respect  to  the  wind.  Thus,  as  the  aircraft  is  flown  in  an  orbit,  • 
the  pilot  only  has  to  vary  the  bank  angle  and  direction  of  flight  of 
the  aircraft.  If  the  aircraft  was  programmed  to  fly  with  a  constant 
ground  track  velocity  in  a  prevailing  wind,  the  pilot  would  have  to 
continually  adjust  the  throttle  in  addition  to  the  bank  angle  and  di¬ 
rection  as  the  aircraft  progressed  around  its  orbit.  Such  a  maneuver 
would  have  che  detrimental  effect  of  increasing  the  yo-yo  motion  of 
the  drogue,  as  noted  in  Case  3  of  Section  A. 2a. 

The  relative  velocity  of  the  aircraft  is  the  vector  summation  of 
the  ground  tr*-  '  'elocity  of  the  aircraft  and  the  negative  of  the  wind 
velocity.  •>  shews  thi3  vector  summation  at  various  points 

along  a  constant  altitude  circular  orbit.  (In  order  to  simplify  ter¬ 
minology,  an  orbit  will  henceforth  be  assumed  to  be  flown  at  a  constant 
altitude  unless  stated  otherwise.)  The  ground  track  velocity  is  tan¬ 
gent  to  the  orbit,  and  the  angle  0  separates  this  velocity  from  the 
positive  Y  axis.  The  angle  \[>  separates  the  aircraft  position  vector 
from  the  positive  X  axis.  The  vector  equations  needed  to  solve  for 
the  ground  track  velocity  are: 

V  -  -  V  i 

w  w 

V  «  -  V  sin8  i+V  cos6  j  . 

t  t  U 
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The  '  jsult..nt  of  these  two  vectors  is: 

V  «  V  +  V  -  (-V  -V„ sine)  3>V„ cos  0  J  . 

ft  W  t  W  t  u 

Taking  the  magnitude  of  each  side  of  the  above  vector  equation  yields 
the  quadratic  equation: 

V*  +  <2Vsine)V,  +  (V  2  -  V  2)  -  0  . 

t  W  t  W  ft 

The  wind  speed,  V  ,  and  the  relative  aircraft  speed,  V  ,  are  specified 

ft 

for  a  particular  flight  test.  The  quadratic  equation  is  solved  for 
the  unknown,  V^.  The  equation  for  V  is  programmed  in  the  computer 

subroutine  that  "flies  the  aircraft".  Once  V  is  calculated,  the  air¬ 
craft  is  flown  at  this  ground  track  velocity  for  an  in  incremental  time 
step.  At  each  new  aircraft  location  V  is  recalculated,  and  in  this 

manner  the  aircraft  proceeds  around  the  orbit. 

If  the  aircraft  flies  in  an  elliptical  orbit,  the  same  equation 

for  V  is  derived.  Figure  3  shows  the  vector  summation  of  V.  and  V  at 
c  t  w 

various  points  along  an  elliptical  orbit.  The  main  difference  between 

the  analysis  of  elliptical  and  circular  orbits  is  that  the  angles  #,  and 

0  are  equal  for  a  circular  orbit  but  not  equal  for  an  elliptical  orbit. 

For  a  circular  orbit  6  *  tan  ^(y/x),  For  an  elliptical  orbit  0  » 

cot  ^(dy/dx). 

A  non-constant  altitude  aircraft  orbit  has  been  proposed  as  a 
method  of  reducing  yo-yo  motion  of  the  drogue  (Crouch,  Baver,  Kahle, 
1970,  p.  36).  When  properly  modified,  a  constant  altitude  circular 
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orbit  transforms  into  a  non-constant  altitude  orbit.  Figure  4  shows 
a  side  view  of  two  non-constant  altitude  orbits.  In  both  orbits,  the 
X  and  Y  coordinates  of  the  aircraft’s  position  are  found  using  the 
computer  subroutine  for  a  constant  altitude  circular  orbit.  The  alti¬ 
tude  of  Orbit  in  Figure  4  is  calculated  from  the  expression: 

Z  -  (X-X  )  tan  +  Zr 
max  C 

This  simple  elliptical  orbit  would  be  impossible  to  fly,  since  the 
vertical  velocity  of  the  aircraft  is  a  very  complicated  function.  How¬ 
ever,  this  particular  orbit  is  a  good  approximation  of  orbits  that  are 
possible  to  fly.  Orbit  2  in  Figure  4  is  ar  orbit  with  a  constant  posi¬ 
tive  vertical  velocity  when  flying  in  the  direction  of  the  prevailing 
wind  and  a  constant  negative  vertical  velocity  when  flying  against  the 
wind.  It  is  more  realistic  to  assume  that  the  orbit  with  constant  ver¬ 
tical  velocities,  i.e.  Orbit  2,  could  be  flown.  Since  both  non-constant 
altitude  orbits  are  very  similar,  the  orbit  that  is  simpler  to  model, 
i.e.  Orbit  1,  is  utilized  in  the  computer  simulation. 

2,3  Drogue  Model 

The  drogue  in  the  computer  simulation  is  modeled  by  a  sphere,  as 
shown  in  Figure  1.  The  final  mass  point  of  the  lumped  mass  model  of 
the  cable  coincides  with  the  center  of  the  sphere.  If  the  cable  is 
assumed  to  attach  to  the  center  of  the  sphere,  all  the  forces  due  to 


the  final  cable  segment  and  due  to  the  sphere  will  act  through  the 
center  of  gravity,  and  the  three  translational  degrees  of  freedom  of 
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the  final  mass  point  will  completely  describe  its  motion.  Thus,  the 
lumped  mass  model  is  modified  at  the  final  mass  point  to  include  the 
weight,  drag,  and  inertia  force  of  the  spherical  drogue  as  forces 
acting  on  the  final  mass  point. 

An  exact  model  of  a  drogue  must  include  rotation.  The  cable  can¬ 
not  be  attached  at  the  center  of  a  sphere  with  a  finite  radius,  so  the 
tension  in  the  cable  produces  a  moment  about  the  center  of  gravity  of 
the  sphere.  For  a  drogue  with  a  more  complex  shape,  the  center  of 
pressure  will  not  coincide  with  the  center  of  gravity,  and  aerodynamic 
forces  will  produce  additional  moments  about  the  center  of  gravity. 
These  moments  produce  rotation,  and  three  additional  equations  of 
motion  are  necessary  to  describe  the  rotation. 

The  equations  of  rotational  motion  could  be  derived,  but  their 
addition  would  not  significantly  improve  the  analysis.  During  SNOWBIRD 
flight  tests,  a  25  lb  conical  drogue  has  been  observed  to  fly  with  its 
axis  always  within  10°  of  the  vercical  in  both  a  zero  wind  environment 
ana  a  prevailing  wind.  A  heavier  drogue  would  have  a  greater  amount  of 
inertia  and  would  align  its  vertical  axis  to  within  a  few  degrees  of 
the  vertical.  Therefore,  the  drogue's  rotational  equations  of  motion 
are  neglected  in  this  analysis. 

2.4  Determination  of  Number  of  Lumped  Masses  in  Model 

If  there  were  an  infinite  number  of  mass  points  in  the  lumped  mass 


model  of  the  trailing  wire  and  if  the  equations  of  motion  could  be 


solved  by  the  computer  with  no  round-off  error,  the  motion  of  the 
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trailing  wire  would  be  described  exactly  by  the  coordinates  of  the 
lumped  masses.  However,  since  only  a  finite  number  of  equations  can 
be  programmed  on  the  computer,  error  in  the  model  is  inevitable.  To 
complicate  the  situation,  the  time  step  used  by  the  Runge-Kutta  inte¬ 
gration  scheme  must  be  very  small  to  keep  the  integration  from  becoming 
unstable.  Therefore,  N,  which  is  tha  number  of  lumped  masses  in  the 
model  of  a  cable,  must  be  minimized  to  prevent  excessive  compute?  usage, 
but  N  must  be  large  enough  to  reduce  the  error  in  the  model  to  an  accep¬ 
table  level. 

Computer  tests  indicate  the  optimum  value  of  N.  The  computer 
program  used  for  the  tests  calculates  the  steady  state  shape  of  the 
trailing  wire  in  a  zero  wind  environment  (Hinnerichs /Crist,  1972). 

During  the  tests,  the  aircraft  flies  in  a  circular  orbit.  Since  there 
is  no  prevailing  wind,  the  cable's  lumped  masses  and  the  drogue  also 
fly  in  circular  orbits.  The  drogue's  orbit  is  specified  and  the  program 
calculates  the  position  of  each  mass  point  leading  to  the  aircraft. 

The  process  cannot  begin  at  the  aircraft,  since  the  tension  in  the  cable 
at  the  aircraft  is  initially  unknown.  As  the  number  of  mass  points  in 
the  lumped-mass  model  increases,  modeling  errors  decrease.  As  N  is  in¬ 
creased,  eventu  *.ly  it  reaches  a  limiting  value,  above  which  a  further 
increase  in  N  results  in  only  a  minute  reduction  in  modeling  errors. 

Both  physical  dimensions  of  the  aircraft  orbit  and  cable  tensions  for  an 
N  greater  than  or  equal  to  this  limiting  N  are  nominal  values  upon  which 
to  base  errors.  An  acceptable  N  produces  errors  no  larger  than  a  few 
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per  cent  of  the  nominal  values. 


Tests  were  performed  for  each  of  the  cable  configurations  analyzed 


in  this  paper.  Table  1  tabulates  modeling  errors  associated  with 


different  values  of  N  for  two  cable  configurations  utilizing  the  SNOW-. 


BIRD  parameters  described  in  Section  3.  The  errors  are  based  on  steady 


state  shape  calculations  for  the  cables  modeled  by  200  lumped  masses. 


For  both  cables  with  N  equal  to  10  the  errors  are  too  large,  and  with  N 


equal  to  20  the  errors  are  acceptable.  As  N  is  increased  beyond  20, 


modeling  errors  are  not  sufficiently  reduced  to  justify  the  increase  in 


computer  time  necessary  to  integrate  the  equations  of  motion.  The 


tensions  are  always  less  for  the  200  limped  mass  model,  so  tensions  in 


the  10  Cases  of  Section  4.2a  are  conservative.  The  progression  angle 


is  defined  to  be  the  angle  between  the  aircraft  and  the  drogue,  with 


its  vertex  at  the  center  of  the  aircraft  orbit.  The  calculations  in 


these  tests  take  only  a  few  seconds  of  central  processor  time  on  the 


Control  Data  Corp  6600  Computer. 


A  time  step  does  not  appear  in  the  steady  state  shape  program,  as 


this  program  does  not  involve  integration.  Thus,  the  optimal  N  should 


be  determined  before  attempting  to  find  the  optimal  time  step. 


2.5  Determination  of  Time  Step 


The  time  step  used  by  the  Runge-Kutta  integration  schen  is  called 


TSTEP.  If  TSTEP  is  too  large,  the  integration  will  become  unstable 


and  the  variab3.es  in  the  equations  of  motion  will  grow  without  bound. 


If  TSTEP  is  too  small,  the  integration  will  be  stable,  but  valuable 
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computer  time  will  be  wasted  and  round-off  errors  could  be  significant. 

The  optimal  TSTEP  is  found  by  choosing  values  for  TSTEP  and  run¬ 
ning  the  computer  simulation.  The  cable  is  initialized  in  its  steady 
state  shape  and  flown  in  a  zero  wind  environment.  The  drogue  and  the 
mess  points  should  not  move  vertically  if  the  integration  is  stable. 

The  largest  TSTEP  for  which  the  integration  is  stable  is  the  best  TSTEP. 
If  TSTEP  is  marginally  stable,  the  integration  may  not  become  unstable 
immediately.  To  insure  that  TSTEP  is  not  marginally  stable,  a  TSTEP 
smaller  than  the  largest  apparently  stable  TSTEP  should  be  chosen.  An 
initial  guess  at  TSTEP  should  be  small  enough  to  integrate  the  highest 
frequency  motion,  which  is  the  longitudinal  vibration  of  the  lumped 
masses.  The  initial  guess  is  given  by: 


TSTEP  - 


10  5  v  2  x  k , 


where:  x  «  period  of  vibration 


m^  *  mass  of  i  point 


*  spring  constant  for  segment  i. 

The  Runge  Kutta  integration  requires  that  TSTEP  “  .016  sec  for 
the  model  of  the  8,000  ft  SNOWBIRD  cable  consisting  of  20  lumped  masses. 
This  is  equivalent  to  .0008  sec/step/mass.  Integrating  with  TSTEP  • 
.0163  takes  approximately  4.8  minutes  of  central  processor  time  on  a 
CDC  6600  Computer  for  each  circular  orbit  flown  in  a  zero  wind  environ¬ 
ment.  Since  the  actual  flying  time  is  about  20.5  sec/orbit,  the  ratio 
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of  central  processor  time  to  aircraft  flight  time  is  approximately 
14:1.  In  a  prevailing  wind,  the  flying  time  will  increase  slightly. 

For  example,  in  a  15  knots  true  air  speed  (KTAS)  wind,  flight  time  is 
24.2  sec/orbit,  and  in  a  30  KTAS  wind,  flight  time  is  30.4  sec/orbit. 

SNOWBIRD  PARAMETERS 

The  parameters  to  be  varied  for  a  complete  parametric  study  of  a 
trailing  wire  system  would  include  the  size  and  shape  of  the  aircraft 
orbit,  the  relative  velocity  of  the  aircraft,  the  physical  properties 
of  the  cable,  and  the  physical  properties  of  the  drogue.  Each  trailing 
wire  configuration  would  be  flown  under  a  variety  of  wind  conditions. 
The  results  of  this  paper  constitute  a  limited  parametric  study  of  the 
SNOWBIRD  system,  a  trailing  wire  system  flight  tested  by  Air  Force 
Academy  personnel. 

3.1  Aircraft  Orbits 

The  aircraft  orbits  used  for  the  computer  simulation  of  the 
SNOWBIRD  system  are  typical  of  the  orbits  flown  by  the  U-10  during 
SNOWBIRD  testing.  The  aircraft  velocity  is  72  knots  true  air  speed 
(KTAS)  relative  to  the  wind  for  both  circular  and  elliptical  orbits. 

For  a  circular  orbit  flown  with  a  zero  wind  condition,  this  aircraft 
velocity  corresponds  to  the  U-10  flying  at  a  30°  batik  angle  with  an 
indicated  speed  of  either  60  knots  indicated  air  speed  (KIAi*>  at  an 
altitude  of  12,825  ft  (4,000  ft  cable)  or  55  KIAS  at  an  altitude  of 
16,900  ft  (measured  from  sea  level)  (8,000  ft  cable).  The  altitude  of 
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the  towing  aircraft  varies  in  order  to  maintain  eimilar  drogue  altitudes 
for  different  length  cables.  The  radius  of  every  circular  orbit  in  the 
simulation  is  800  ft.  The  orbital  parameters  for  a  circular  orbit 
flown  with  a  zero  wind  condition  are  calculated  from  the  equations: 


where: 


R  -  T* 

* 

$  *  angular  velocity  of  aircraft 


V  ■  tangential  or  ground  track  velocity  of  aircraft 
X  »  bank  angle  of  aircraft 


R  =>  radius  of  circular  orbit 


g  “  acceleration  of  gravity 

(Liller,  1969) 

If  a  prevailing  wind  exists,  X  must  continually  be  corrected  in  order 
to  fly  a  circular  orbit  and  both  V  and  d>  will  vary;  but  the  relative 

aircraft  velocity,  the  orbit  altitude,  and  R  are  the  same  as  in  a  zero 
wind  environment.  All  elliptical  orbits  in  the  simulation  have  a  semi¬ 
major  axis  of  1200  ft  and  a  semi -minor  axis  of  800  ft. 

The  circular  orbit  in  Figure  2  and  the  elliptical  orbit  in  Figure 
3  represent  SNOWBIRD  orbits. 

3.2  Cable  Properties 

The  following  cable  properties  were  measured  for  the  SNOWBIRD 


cable: 


diameter  «  .048  in. 
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Vortex  shedding  has  also  been  omitted  from  most  studies  of  cable 
systems  under  hydrodynamic  loading  (Casarella  and  Parsons,  1970,  p.  35) 
3.3  Drogue  Properties 

The  following  properties  were  assumed  for  the  spherical  drogue: 
Radius  **  .42  ft. 

Weight  «•  25  lb. 

In  one  simulated  flight  test  the  drogue  weight  was  increased  to  100  lb 
while  the  radius  remained  constant. 


Aerodynamic  lift  and  drag  coefficients  for  a  spherical  drogue  are 


expressed  in  the  form:  C  *  0. 

Li 


.6  . 


Lift  and  drag  forces  are  given  by:  L  ■  0. 

D  »  CDqA 

2 

where:  A  ■  cross  sectional  area  of  sphere  (ir  x  radius  ). 

The  air  density  for  all  aerodynamic  force  calculations  is  taken  from  a 
curve  fit  of  the  NACA  Standard  Atmosphere  Tables. 

3.4  Wind  Conditions 

Two  types  of  wind  profiles  were  used  in  the  computer  analysis. 

The  simplest  and  most  frequently  used  is  the  constant  wind  profile. 
This  profile  is  used  when  comparing  dissimilar  cable  configurations. 
The  constant  wind  profile  assures  the  analyst  that  variations  in  the 
cable  shapes  and  motions  are  not  due  to  a  change  in  shape  of  the  wind 
profile,  but  instead  are  a  consequence  of  differences  in  the  physical 
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parameters  of  the  cable  configurations.  The  power  law  profile  is  used 
to  approximate  the  wind  profile  at  low  altitudes.  The  power  law  is 
stated  as: 

r 

h. 


m 


v  +  V 
V2  V1 


where:  ■  wind  speed  at  altitude  h ^ 

»  wind  speed  at  altitude  h^ 
m  -  .35 

(Liller,  1969,  p.  38) 

Since  SNOWBIRD  test  flights  were  flown  at  a  location  where  ground  level 
was  8,900  ft,  in  the  computer  simulation  h^  *  12,900  ft  and  h£  ranges 

from  8,900  ft  to  12,900  ft.  At  altitudes  greater  than  h^,  the  wind 

profile  becomes  a  constant  profile. 

The  wind  speed  for  the  computer  simulation  took  on  values  of  0, 

15,  and  30  KTAS . 

DISCUSSION  OF  RESULTS 

4.1  Determination  of  Cable  Lengths 

Figure  6  is  a  plot  of  the  radius  of  the  drogue  orbit  for  various 

cable  lengths  flown  in  a  zero  wind  environment.  The  aircraft  orbit  is 

the  same  for  each  cable  length,  and  is  the  same  size  as  the  circular 

orbit  in  Figure  2.  Since  V  is  0,  both  V  and  V,,,  are  72  KTAS.  The 

w  a  T 

drogue  altitude  is  9,100  ft.  The  aircraft  altitude  ranges  from  18,940 
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ft  for  a  10,000  ft  cable  to  9,420  ft  for  a  500  ft  cable. 

In  order  to  study  the  effect  that  the  cable  length  has  upon  the 
motion  of  a  trailing  wire  system  in  a  prevailing  wind,  two  lengths  of 
cable  from  Figure  6  were  chosen  to  be  flight  tested.  A  4,000  ft  cable 
was  chosen,  since  for  a  shorter  cable  the  drogue  orbital  radius  is  in¬ 
tolerably  large.  An  8,000  ft  cable  was  chosen,  since  increasing  the 
cable  length  beyond  8,000  ft  does  not  significantly  reduce  the  drogue 
orbital  radius. 

4.2  Comparison  of  Flight  Configurations 

A  computer  simulation  consisting  of  twelve  different  orbiting 
trailing  wire  configurations  has  been  arranged  into  ten  representative 
cases  (Case  1-10)  in  a  movie.  (Trailing  Wire  Antenna  Configurations 
Using  Snowbird  Cable  Towed  by  U-10  Aircraft)  Individual  frames  of  the 
movie  are  presented  in  this  report.  Time  exposures  of  each  case  follow 
in  the  analysis. 

4.2a  Individual  Frames 


Five  individual  frames  of  Cases  1-10  appear  in  Figures  7a  through 
16e.  For  each  case  the  sequence  of  frames  labeled  ,!a"  through  "e" 
presents  similar  information.  Figures  7-16a  are  title  frames  describing 
the  case  and  providing  a  legend  in  order  to  distinguish  between  differ¬ 
ent  cables.  Figure  7-16b  show  three  orthogonal  projections  of  the 
cables.  The  aircraft  is  located  at  the  point  on  the  orbit  coincident 
with  the  positive  X  axis.  The  cables  appear  to  be  smooth  curves,  but 
they  actually  consist  of  20  straight  line  segments.  Figures  7-16c  are 
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plan  views,  i.e.,  projections  onto  a  horizontal  plane,  of  the  aircraft 
and  drogue  orbits.  The  aircraft  and  drcgue  fly  counterclockwise  in 
their  orbits.  Figure  7-16d  are  projections  of  the  drogue  orbit  onto  a 
vertical  plane  containing  the  X  axis.  These  figures  depict  the  verti¬ 
cal  oscillation  of  the  drogue,  or  the  drogue  yo-yo.  Figures  7-16e  plot 
the  tensions  in  the  cable  segment  attached  to  the  aircraft.  This  seg¬ 
ment  is  particularly  interesting,  since  cable  tensions  are  generally 
greatest  at  the  aircraft.  In  Figures  7-16d  and  Figures  7-16e,  the  yo¬ 
yo  and  the  tension  are  plotted  for  four  identical  orbits  of  the  air¬ 
craft.  The  motion  of  the  cable  is  identical  for  each  cycle  because  the 
cable  has  been  flying  long  enough  to  attain  a  steady  state  shape. 

Table  2  depicts  a  tabular  summary  of  the  information  contained  in 
Figures  7a  through  16e. 

Two  suspended  cables  in  the  X-Z  projections  of  Figures  9-12b  pre¬ 
sent  additional  information  for  Cases  9-12.  These  suspended  cables, 
being  in  equilibrium,  are  motionless.  One  cable  is  suspended  from  the 
center  of  the  aircraft  orbit,  while  the  second  cable  hangs  from  the 
position  on  the  orbit  coincident  with  the  positive  X  axis.  The  en¬ 
velope  for  the  motion  of  the  towed  cable  consists  of  the  area  between 
the  two  suspended  cables  in  the  X-Z  projection.  The  entire  cable  will 
not  remain  within  the  envelope  during  a  complete  aircraft  orbit.  Be¬ 
cause  the  wind  blows  in  the  positive  x  direction,  the  uppermost  section 
of  the  cable  is  located  outside  the  envelope  when  the  aircraft  has  a 
negative  X  coordinate,  while  the  lowermost  section  of  cable  never  leaves 


the  envelope. 


Case 

Time (Sec) 
Cycle 

Drogue 

Drogue  Average 

Peak  to  Peak  Ascending 
Yo-Yo  (Ft)  Speed 

(ft/sec) 

Drogue 

Average 

Descending 

Speed 

(ft/sec) 

Peak  Tension 
in  Cable 
at  Aircraft 

1 

(Blue) 

41.0 

0. 

0. 

0. 

63.8 

1 

(Black) 

41.0 

0. 

0. 

0. 

43.9 

2 

(Blue) 

52.2 

129. 

9.0 

10.8 

69.9 

2 

(Black) 

52.2 

179. 

13.2 

14.5 

48.7 

3 

47.9 

1,247. 

38.9 

78.7 

144.9 

4 

42.9 

682. 

26.3 

40.1 

89.1 

5 

53.4 

736. 

20.2 

42.3 

92.6 

6 

54.1 

990. 

32.0 

42.7 

78.8 

7 

(Blue) 

Same  as 

Case  4 

7 

(Black) 

42.9 

608. 

24.2 

34.2 

58.2 

8 

(Blue) 

42.9 

332. 

13.0 

19.1 

163.1 

8 

(Black) 

Same  as 

Case  4 

9 

(Blue) 

42.9 

338. 

15.2 

22.3 

76.3 

9 

(Black) 

Same  as 

Case  4 

10 

(Blue) 

42.9 

603. 

23.3 

36.4 

86.2 

10 

(Black) 

Same  as 

Case  4 

TABLE  2 


SUMMARY  OF  INFORMATION  CONTAINED 
IN  FIGURES  7a-16e 
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This  case  illustrates  the  advantages  of  towing  an  8,000  ft  cable 
compared  to  a  4,000  ft  cable  when  flying  a  circular  orbit  in  a  zero 
wind  environment. 
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Figure  7b 

The  projections  show  that  the  8,000  ft  cable  has  much  longer  hang¬ 
ing  chan  the  4,000  ft  cable,  i.e.,  a  straight  line  drawn  from  the  air¬ 
craft  to  the  drogue  provides  a  better  approximation  to  the  shape  of  the 
4,000  ft  cable.  The  progression  angles  are  290°  for  the  8,000  ft  cable 
and  195°  for  the  4,000  ft  cable.  Thus  the  8,000  ft  cable  has  more  coil¬ 
ing  and  therefore  more  slack  than  the  4,000  ft  cable. 

Figure  7c 

The  plan  view  of  the  drogue  orbit  for  both  cables  is  a  circle.  If 

the  aircraft  flies  in  a  circular  orbit  with  no  prevailing  wind,  every 

point  on  the  cable  will  travel  in  a  circular  orbit. 

Since  the  aircraft  and  the  drogue  both  have  the  same  angular  velo¬ 
city,  the  tangential  velocity  of  the  drogue  wil3  be  less  than  the  air¬ 
craft  tangential  velocity.  The  tangential  velocities  for  this  case  are: 
VT  aircraft  =  122.  ft/sec 

drogue  for  8,000  ft  cable  =>  21.  ft/sec 

VT  drogue  for  4,000  ft  cable  »  46.  ft/sec 


Figure  7d 

For  an  aircraft  flying  in  a  circular  orbit  in  the  absence  of  a 
prevailing  wind,  each  point  on  the  cable  orbits  at  a  constant  altitude. 
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Therefore  the  drogue  yo-yo  is  zero.  (In  this  context  the  term  yo-yo 
applies  to  the  peak  to  peak  displacement  of  the  drogue's  vertical 
oscillation.) 

Figure  7e 

The  peak  cable  tensions  are  a  few  pounds  less  than  the  system  dead 
weights  for  both  the  8,000  ft  and  4,000  ft  cables.  The  dead  weight  of 
a  trailing  wire  system  consists  of  the  weight  of  the  cable  and  the 
weight  of  the  drogue.  The  dead  weight  of  the  SNOWBIRD  cable  wich  a 
25  lb  drogue  is  65  lb  for  the  8,000  ft  cable,  and  45  lb  for  the  4,000 
ft  cable. 

Both  the  8,000  ft  and  4,000  ft  systems  in  CASE  1  provide  a  plat¬ 
form  that  could  be  used  for  many  practical  applications.  Tensions  are 
not  near  the  300  lb  breaking  strength  of  the  particular  cable  used. 

Both  the  tangential  velocity  and  the  drogue  plan  area  are  less  for  the 
8,000  ft  cable,  so  it  would  be  easier  to  utilize. 

Flight  conditions  have  been  idealized  for  this  case,  since  there  is 
no  wind  and  the  aircraft  flies  in  a  circular  orbit.  CASE  2  studies  the 
same  two  cable  systems;  however,  the  aircraft  flies  in  an  elliptical 
orbit. 

CASE  2 

This  case  illustrates  the  advantages  of  towing  an  8,000  ft  cable 
compared  to  a  4,000  ft  cable  when  flying  an  elliptical  orbit  in  a  zero 
wind  environment. 
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Figure  8b 

The  projected  views  of  this  case  are  very  similar  to  the  projected 
views  of  CASE  1  (Figure  7b).  The  alack  is  again  greater  for  the  8,000 
ft  cable. 

Figure  8c 

The  plan  view  of  the  drogue  orbit  for  both  cable  systems  resembles 
an  ellipse.  The  plan  views  of  both  the  aircraft  and  drogue  elliptical 
orbits  circumscribe  the  corresponding  plan  views  of  the  aircraft  and 
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drogue  circular  orbits  of  CASE  1  (Figure  ?c).  However,  if  the  aircraft 
circular  orbit  circumscribed  the  aircraft  elliptical  orbit,  the  drogue 
circular  orbit  would  circumscribe  the  drogue  elliptical  orbit. 

Figure  8d 

Since  the  aircraft  orbit  is  not  circular,  the  drogue  oscillates 
vertically.  The  two  crests  in  a  complete  cycle  of  the  drogue's  motion 
occur  as  the  aircraft  crosses  the  Y  axis,  and  the  two  troughs  per  cycle 
occur  as  the  aircraft  crosses  the  X  axis.  As  the  aircraft  flies  coun¬ 
terclockwise  around  the  orbit  from  the  positive  X  axis  to  .the  positive 
Y  axis,  the  slack  in  the  cable  decreases  and  the  drogue  rises.  The 
more  slack  that  exists  in  the  system,  the  less  the  drogue  will  rise. 
Since  there  is  more  alack  in  the  8,000  ft  cable  than  in  the  4,000  ft 
cable,  the  yo-yo  is  less  for  the  8,000  ft  cable. 

Figure  8e 

The  peak  tensions  for  both  cable  systems  are  less  than  10%  greater 
than  their  respective  dead  weights.  Flying  the  system  in  a  prevailing 


itw^W’wJw 


wind  will  produce  much  greater  tensions. 

The  8,000  ft  cable  has  better  flight  characteristics  than  the 
4,000  ft  cable  when  flying  an  elliptical  orbit  in  a  zero  wind  environ¬ 
ment.  Both  the  horizontal  and  vertical  displacements  and  velocities 
are  less  for  the  8,000  ft  cable.  Tensions  in  this  case  are  so  low  that 
they  are  not  critical  in  the  design  of  either  system. 

Comparing  CASE  1  with  CASE  2,  one  concludes  that  in  a  zero  wind 
condition,  the  aircraft  should  be  flown  in  a  circular  orbit  to  obtain 
optimal  flight  characteristics  of  the  drogue.  Because  the  aircraft 
cannot  be  flown  in  a  perfect  orbit,  the  length  of  the  cable  is  important. 
As  the  length  of  the  cable  increases,  the  yo-yo  motion  of  the  drogue 
decreases.  This  has  been  observed  in  SNOWBIRD  flight  tests. 

CASE  3 


This  case  demonstrates  the  detrimental  effect  that  a  30  KTAS  pre¬ 
vailing  wind  has  upon  the  yo-yo  of  the  drogue  and  the  cable  tensions  for 
an  8,000  ft  cable. 

Figure  9b 

In  the  X-Z  projection  the  towed  cable  remains  entirely  within  the  * 
envelope  formed  by  the  two  suspended  cables.  The  trailing  2,000  ft  of  l 
the  cable  hangs  nearly  vertical.  * 

Figure  9c 

The  plan  view  of  the  drogue  orbit  has  a  smaller  area  in  this  case 
than  in  CASE  1  (Figure  7c),  where  the  aircraft  flies  in  a  similar  cir¬ 
cular  orbit  but  with  no  prevailing  wind.  The  drogue  orbits  downwind  of 
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the  aircraft  at  a  distance  slightly  greater  than  one-half  the  cable 
length.  Because  the  aircraft  flies  counterclockwise  around  the  orbit 
and  spends  more  time  flying  upwind  than  downwind,  the  drogue  orbits  on 
the  positive  Y  side  of  the  X  axis. 

Figure  9d 

This  much  yo-yo  could  be  very  detrimental  to  the  effectiveness  of 
the  trailing  wire  system.  The  average  drogue  velocity  during  descent, 
78.7  ft/sec,  makes  it  very  difficult  to  land  the  drogue  softl-  in  a  30 
KTAS  wind. 

The  highest  altitude  of  the  drogue  occurs  when  the  aircraft  is 
furthest  upwind,  and  its  lowest  altitude  occurs  when  the  aircraft  is 
furthest  downwind. 

The  yo-yo  for  the  same  system  flown  under  the  same  wind  condition 
but  with  Vfc  constant  and  equal  to  72  KTAS  is  1,300  ft. 

Figure  9e 

The  peak  tension  is  123%  greater  than  the  dead  weight  of  the  system, 
and  very  near  the  300  lb  breaking  strength  of  the  SNOWBIRD  cable. 

As  the  aircraft  flies  upwind,  the  tension  in  the  cable  continually 
decreases.  When  the  aircraft  flies  downwind,  it  doubles  back  on  the 
cable  and  the  drogue  falls.  Tensions  in  the  cable  decrease  even  further 
as  the  drogue  falls.  If  the  drogue  were  free-falling,  its  velocity 
would  increase  until  it  attained  a  terminal  velocity  of  289  ft/sec. 
Restraining  tension  in  the  cable  causes  the  maximum  drogue  vertical 
velocity  to  be  less  than  100  ft/sec.  Just  prior  to  the  aircraft 
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crossing  the  positive  X  axis,  the  drogue  has  pulled  the  cable  taut,  and 
a  spike  in  the  cable  tension  decelerates  the  drogue.  A  stress  wave 
travels  up  the  cable  from  the  drogue  to  the  aircraft.  The  velocity  of 
the  stress  wave  is  */f7p  or  10,902  ft/sec.  It  takes  0.734  sec  for  the. 
stress  wave  to  travel  up  8,000  ft  of  cable  to  the  aircraft. 

The  yo-yo  of  the  trailing  wire  system  with  a  25  lb  drogue  in  a  30 
KTAS  wind  would  prohibit  its  use  in  many  applications.  Tensions  are  so 
very  near  the  breaking  strength  of  the  SNOWBIRD  cable  that  additional 
gust  loading  at  the  drogue  could  conceivably  result  in  failure  of  the 
cable. 

The  30  KTAS  wind  at  a  constant  profile  creates  an  extremely  poor 
environment  in  which  to  operate  the  trailing  wire  system.  Normally 
wind  speeds  will  not  be  so  high,  since  the  system  is  flown  near  the 
ground.  The  power  law  profile  (explained  in  the  "Wind  Conditions" 
section)  used  in  CASE  10  gives  a  better  approximation  to  the  actual 
wind  profile  at  altitudes  within  4,000  ft  of  the  ground,  and  the  trail¬ 
ing  wire  system  behaves  better  with  a  power  law  wind  profile. 

CASES  4-10  study  methods  used  to  reduce  yo-yo  and  peak  tensions 
when  flying  in  a  prevailing  wind.  The  maximum  wind  speed  for  these 
cases  is  15  KTAS. 

CASE  4 

Thi3  case  demonstrates  the  detrimental  effect  that  a  15  KTAS  pre¬ 
vailing  wind  has  upon  the  yo-yo  of  the  drogue  and  the  cable  tensions 
for  an  8,000  ft  cable. 
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In  the  X-Z  projection  the  towed  cable  remains  entirely  within  the 
envelope  formed  by  the  two  suspended  cables.  The  trailing  4,000  ft  of 
the  cable  hangs  nearly  vertical. 

Figure  10c 

The  plan  view  of  the  drogue  orbit  is  even  smaller  in  a  15  KTAS  wind 
than  in  the  30  KTAS  wind  (Figure  9c),  and  for  all  practical  purposes, 
the  motion  of  the  drogue  is  entirely  vertical. 

Figure  lOd 

Both  the  average  descending  velocity  and  yo-yo  of  the  drogue  are 
approximately  50%  of  the  corresponding  values  for  a  30  KTAS  wind  (Figure 
9d)  .  It  is  still  extremely  difficult  to  land  the  drogue  without  its 
crashing . 

Figure  lOe 

The  peak  tension  is  only  37%  greater  than  the  dead  weight  of  the 
system,  or  less  than  one-half  of  the  breaking  strength  of  the  cable. 

The  tensions  are  less  than  in  a  30  KTAS  wind  (Figure  9e) ,  since  the 
drogue  descending  velocity  is  less. 

The  average  descending  velocity  and  yo-yo  of  the  drogue  and  the 
peak  tension  in  the  cable  caused  by  flying  in  a  15  KTAS  wind  are  all 
approximately  50%  of  the  corresponding  values  for  a  30  KTAS  wind. 

CASE  5 


This  case  studies  an  8,000  ft  cable  towed  in  an  elliptical  orbit 
with  the  major  axis  of  the  ellipse  perpendicular  to  the  direction  of 
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a  15  KTAS  prevailing  wind. 

Figure  lib 

The  projections  are  nearly  identical  to  the  projections  in  CASE  4, 
where  a  circular  orbit  is  flown  under  the  same  wind  conditions  (Figure 
10b).  In  the  X-Z  projection  the  towed  cable  is  entirely  within  the 
envelope  formed  by  the  two  suspended  cables. 

Figure  11c 

The  plan  view  of  the  drogue  orbit  centers  about  the  same  location 
as  the  plan  view  in  CASE  4  (Figure  10c) .  The  area  of  the  plan  view  here 
is  only  slightly  greater  than  the  area  in  CASE  4. 

Figure  lid 

The  yo-yo  is  7.9%  greater  than  the  yo-yo  in  Ca.SE  4  (Figure  lOd)  . 
However,  since  the  path  of  the  elliptical  orbit  is  longer  than  the  path 
of  the  circular  orbit,  the  drogue  average  descending  speed  is  only  5.5% 
greater  in  this  case  than  in  CASE  4.  The  drogue  average  ascending  speed 
is  less  in  this  case  than  i!  CASE  4. 

There  are  four  distinct  cycles  of  yo-yo  motion.  In  CASE  2,  where 
an  identically  shaped  elliptical  orbit  is  flown  in  a  zero  wind  environ¬ 
ment,  there  are  eight  cycles  of  yo-yo  motion.  Thus,  the  yo-yo  for  this 
case  is  primarily  due  not  to  the  shape  of  the  oibit,  but  to  the  le*  ;th 
of  time  spent  flying  into  the  prevailing  wind. 

Figure  lie 

The  peak  tension  is  only  2.5  lb  greater  than  the  peak  tension  in 
CASE  4  (Figure  lOe) . 
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The  drogue  flight  characteristics  and  the  cable  tensions  for  this 
case  are  only  slightly  less  desirable  than  those  for  CASE  4.  From  this 
case  one  concludes  that  irregularities  in  the  shape  of  an  orbit  causing 
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it  to  be  noncircular  will  not  limit  the  application  of  the  trailing 
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wire  system,  if  the  irregularities  only  increase  the.  dimension  of  the 
orbit  perpendicular  to  the  prevailing  wind  by  one-half  the  diameter  of 
the  desired  circular  orbit. 
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CASE  6 
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This  case  studies  an  8,000  ft  cable  towed  in  an  elliptical  orbit 
with  the  major  axis  of  the  ellipse  parallel  to  the  direction  of  a  15 
KTAS  prevailing  wind. 

Figure  12b 

Again  the  projections  are  nearly  identical  to  the  projections  in 
CASE  4  (Figure  10b),  and  in  the  X-Z  projection  the  towed  cable  is  en¬ 
tirely  within  the  envelope  formed  by  the  two  suspended  cables. 

Figure  12c 

The  plan  view  o£  the  drogue  orbit  is  approximately  200  ft  further 
downwind  than  the  orbit  in  CASE  4  (Figure  10c) .  The  area  of  the  plan 
view  is  only  slightly  greater  than  the  area  in  CASE  4, 

Figure  12d 

The  yo-yo  is  45.2%  greater  than  the  yo-yo  in  CASE  4  (Figure  lOd) . 
Since  the  path  of  the  elliptical  orbit  is  longer  than  the  path  of  the 
circular  orbit,  the  drogue  average  descending  speed  is  only  6.5%  greater 
in  this  case  than  in  CASE  4. 
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As  in  CASE  5,  there  are  four  distinct  cycles  of  yo-yo  motion,  and 
the  yo-yo  is  primarily  due  not  to  the  shape  of  the  orbit,  but  to  the 
length  of  time  spent  flying  into  the  prevailing  wind. 

Figure  12e 

The  peak  tension  is  10.3  lb  less  than  the  peak  tension  in  CASE  A 
(Figure  lOe) . 

The  only  factor  making  this  case  less  desirable  than  CASE  A  is  the 
increase  in  yo-yo.  Since  the  drogue  average  descending  spe^d  is  com¬ 
parable  to  the  speed  in  CASE  A,  the  drogue  in  this  case  would  not  be  any 
more  difficult  to  land  than  the  drogue  in  CASE  A.  The  decrease  in  the 
peak  tension  in  the  cable  is,  of  course,  desirable. 

From  this  case  one  concludes  that  irregularities  in  the  shape  of 
an  orbit  causing  it  to  be  noncircular,  would  limit  the  application  of  a 
trailing  wire  system  only  if  yo-yo  had  to  be  kept  to  an  absolute  mini¬ 
mum.  Significant  irregularities  are  those  that  cause  the  dimension  of 
the  orbit  parallel  to  the  prevailing  wind  to  be  increased  by  one-half 
the  diameter  of  the  desired  circular  orbit. 

CASES  A,  5  and  6  together  infer  that  in  a  prevailing  wind:  1)  a 
circular  orbit  minimizes  the  yo-yo,  average  descending  vertical  velocity, 
and  orbital  plan  area  of  the  drogue;  and  2)  irregularities  in  the  shape 
of  the  orbit,  such  that  no  dimension  of  the  orbit  is  greater  than  one 
and  one-half  times  the  diameter  of  the  desired  circular  orbit,  will  not 
limit  i.  application  of  a  trailing  wire  system.  The  analysis  of  time 
exposures  of  the  cable  motion  will  further  conclude  that  the  smallest 
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circular  orbit  is  the  most  desirable. 

CASE  7 

This  case  compares  an  8,000  ft  cable  with  a  4,000  ft  cable,  when 
towed  in  identical  circular  orbits  in  a  15  KTAS  prevailing  wind. 

Figure  13b 

The  cable  projections  are  significantly  different  from  those  of 
CASE  1  (Figure  7b) ,  in  which  the  same  cables  are  flown  in  a  zero  wind 
environment . 

Figure  13c 

The  plan  areas  of  the  drogue  orbits  in  this  case  are  less  than  the 
areas  in  CASE  1  (Figure  7c).  The  area  for  the  8,000  ft  cable  is  sig¬ 
nificantly  less  than  for  the  4,000  ft  cable  in  both  cases.  Therefore, 
horizontal  drogue  velocities  are  significantly  less  for  the  8,000  ft 
cable. 

Figure  13d 

The  yo-yo  for  the  4,000  ft  cable  is  10.8%  less  than  the  yo-yo  for 
the  8,000  ft  cable.  The  drogue  average  descending  velocity  is  14.7% 
less  for  the  4,000  ft  cable,  so  the  drogue  is  slightly  easier  to  land 
when  flown  on  a  4,000  ft  cable. 

Figure  13e 

For  the  8,000  ft  cable  the  peak  tension  is  37%  greater  than  its 
system  dead  weight,  and  for  the  4,000  ft  cable  the  peak  tension  is  29.4% 
greater  than  its  deadweight.  The  tension  increases  more  for  the  8,000 
ft  cable,  since  its  drogue  average  descending  velocity  is  greater. 
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The  yo-yo  and  average  descending  velocity  of  the  drogue  and  the 
peak  cable  tension  are  all  less  for  the  4,000  ft  cable,  but  the  hori¬ 
zontal  drogue  motion  is  less  for  the  8,000  ft  cable. 

CASE  8 

This  case  compares  two  trailing  wire  systems,  both  utilizing  8,000 
ft  cables  and  towed  in  identical  circular  orbits  in  a  15  KTAS  prevailing 
wind.  One  system  has  a  25  lb  drogue  and  the  other  a  100  lb  drogue. 
Although  the  drogues  are  of  different  weights,  they  have  the  same  dimen¬ 
sions  . 

Figure  14b 

The  cable  with  the  100  lb  drogue  hangs  almost  directly  beneath  the 
aircraft. 

Figure  14c 

The  plan  area  of  the  drogue  orbit  is  much  larger  for  the  cable 
with  the  100  lb  drogue.  If  the  two  systems  were  flown  in  still  air, 
the  drogue  orbital  radius  would  also  be  larger  for  the  system  with  the 
heavier  drogue.  The  plan  area  and  location  of  the  100  lb  drogue's  orbit 
resemble  the  plan  area  and  location  of  the  4,000  ft  cable  with  a  25  lb 
drogue  in  CASE  7  (Figure  13c).  Therefore,  both  shortening  the  cable 
and  increasing  the  weight  of  the  drogue  increase  the  size  of  the  plan 
view  of  the  drogue  orbit. 

Figure  14d 

The  yo-yo  of  the  100  lb  drogue  is  51.3%  less  than  the  yo-yo  of  the 
25  lb  drogue,  and  its  average  descending  velocity  is  52.4%  less  than 
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for  the  25  lb  drogue. 

Figure  14e 

The  peak  tension  in  the  cable  with  the  100  lb  drogue  is  23  lb 
greater  than  its  system  dead  weight,  while  in  the  cable  with  the  25 
lb  drogue,  the  peak  tension  is  24.1  lb  greater  than  its  system  dead 
weight.  The  increase  in  tension  above  the  dead  weight  is  nearly  the 
same,  because  for  the  heavier  system,  the  smaller  drogue  descending 
velocities  are  counteracted  by  the  larger  mass  of  the  drogue.  The  peak 
tension  for  the  system  incorporating  the  100  lb  drogue  is  dangerously 
near  the  cable's  300  lb  breaking  point. 

The  yo-yo  and  average  descending  velocity  cf  the  drogue  are  less 
for  the  system  with  the  100  lb  drogue.  The  peak  cable  tension  and  hori¬ 
zontal  drogue  motion  are  less  for  the  system  with  the  25  lb  drogue, 

CASE  9 

This  case  compares  two  identical  8,000  ft  trailing  wire  systems 
flown  in  a  15  KTAS  wind.  One  is  towed  in  a  constant  altitude  circular 
orbit,  and  the  other  in  the  nonconstant  altitude  orbit  labeled  "Orbit  1" 
in  Figure  4. 

Figure  15b 

The  projections  of  the  two  systems  are  very  similar. 

Figure  15c 

The  plan  views  of  the  drogue  orbits  are  also  nearly  identical. 

Figure  15d 

The  yo-yo  and  drogue  average  descending  velocity  of  the  system 


employing  a  nonconstant  altitude  orbit  are  43.1%  and  44.4%  less  than  the 
corresponding  values  for  the  system  with  a  constant  altitude  orbit. 

The  difference  in  altitude  between  the  highest  and  lowest  points 
in  the  nonconstant  altitude  orbit  amounts  to  300  ft.  The  yo-yo  for  this 
orbit  is  294  ft  less  than  for  the  constant  altitude  orbit.  Analysis  of 
Section  4.2b  "Time  Exposures"  will  explain  why  this  is  such  an  excellent 
method  of  reducing  yo-yo. 

Figure  15e 

The  peak  cable  tension  is  14.4%  less  for  the  nonconstant  altitude 
orbit,  because  the  drogue  descending  velocities  are  less. 

The  yo-yo  and  average  descending  velocity  of  the  drogue  and  the 
peak  cable  tension  are  all  less  for  the  system  flown  ir.  a  nonconstant 
altitude  orbit.  The  horizontal  motion  remains  approximately  the  same 
for  both  systems.  Of  the  methods  studied  in  this  report,  flying  in  a 
nonconstant  altitude  aircraft  orbit  is  optimal. 

CASE  10 


This  case  compares  two  8,000  ft  trailing  wire  systems  flown  in 
identical  circular  orbits  with  a  15  KTAS  wind  at  the  aircraft's  orbit. 
For  one  system  the  wind  profile  is  constant  and  for  the  other  a  power 
law  wind  profile  is  in  effect. 

Figure  16b 

The  cable  loaded  by  the  wind  in  a  power  law  profile  is  exposed  to 
lower  wind  velocities  than  the  cable  loaded  by  the  wind  in  a  constant 
profile  as  explained  in  Section  3.4;  hence,  it  is  not  blown  as  far  down¬ 
wind. 
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Figure  16c 

The  reduced  wind  velocities  result  in  a  larger  plan  area  of  the 
drogue  orbit. 

Figure  16d 

The  yo-yo  and  drogue  average  descending  velocity  are  less  by  11,6% 
and  9.3%  respectively  for  the  system  incorporating  the  power  law  wind 
profile. 

Figure  16e 

The  peak  cable  tension  is  reduced  by  3.3%  when  the  power  law  wind 
profile  is  assumed. 

It  has  been  reported  that  the  wind  acting  on  the  lower  section  of 
the  cable  influences  the  motion  of  the  drogue  to  a  greater  extent  than 
the  wind  on  the  upper  section  of  the  cable.  (Liller,  1969,  p.  59)  The 
reduced  wind  velocities  attributed  to  the  power  law  profile  occur  near 
the  ground  and  affect  the  lower  section  of  the  cable.  When  the  power 
law  wind  profile  is  used,  the  vertical  motion  is  reduced  by  roughly  10% 
in  a  15  KTAS  wind.  Assuming  that  the  reduction  remains  10%  for  a  30 
KTAS  wind,  the  reduced  vertical  motion  derived  from  Case  3  would  still 
be  very  undesirable. 

4 . 2b  Time  Exposures 

Time  exposures  of  Cases  1-9  appear  in  Figures  18-26.  The  exposures 
show  the  three  orthogonal  projections  that  appear  in  the  computer  simu¬ 
lated  movie. 
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CASE  1 


Figure  18  depicts  surfaces  swept  out  by  the  cables  when  there  is 
no  wind.  The  funnel  shaped  surfaces  are  symmetric  about  the  X  axis. 

Both  surfaces  attain  their  smallest  cross-sectional  diameter  above  the 
lower  end  of  the  cable  and  flare  outward  beneath  this  point  due  to  the 
centrifugal  force  of  the  drogue. 

CASE  2 

The  surfaces  in  Figure  19  are  also  funnel-sheped.  In  this  case 
cross  sections  of  the  surfaces  in  the  X-Y  plane  are  ellipses  and  not 
circles  as  in  Case  1. 

CASE  3 

The  X-Z  projection  of  Figure  20  shows  that  in  a  30  KTAS  wind  the 
lowest  5,500  ft  of  cable  remain  within  the  envelope  formed  by  the  sus¬ 
pended  cables  during  an  entire  orbit.  There  is  no  flare  at  the  bottom 
of  the  surface.  Thus  the  plan  area  of  the  drogue  orbit  is  smaller  if 
the  trailing  wire  system  is  flown  in  a  prevailing  wind,  than  if  it  is 
flown  in  a  zero  wind  condition. 

Examination  of  the  towed  cable's  motion  with  respect  to  the  envelope 
formed  by  the  suspended  cables  suggests  a  method  of  calculating  the  yo-yo 
and  the  average  descending  and  ascending  speeds  of  the  drogue  without 
having  to  integrate  the  dynamic  equations  of  motion.  The  lowest  section 
of  cable  moves  in  a  path  almost  parallel  to  the  suspended  cables.  When 
the  aircraft  flies  upwind,  the  cable  travels  up  the  path,  and  as  the 
aircraft  flies  downwind,  the  cable  travels  back  down  a  nearly  identical 
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path.  This  phenomenon  occurs  because  the  lumped  mass  model  of  the  cable 
is  perfectly  flexible,  i.e.,  the  cable  can  support  no  bending  moment. 

In  actuaxity  the  cable  can  support  only  a  negligible  bending  moment,  so 
the  model  is  a  very  good  approximation.  The  uppermost  section  of  the 
cable  sweeps  out  a  funnel  shaped  surface  as  the  plane  orbits,  thus  this 
section  of  the  cable  does  not  remain  parallel  to  the  suspended  cables. 

If  the  junction  between  the  uppermost  and  lowermost  sections  of  the 
cable  was  located  on  the  Z  axis,  there  would  be  no  yo-yo  motion  of  the 
drogue.  Since  the  junction  is  downwind  of  the  Z  a.ris,  the  distance  from 
the  aircraft  to  the  junction  is  greater  when  the  aircraft  has  a  negative 
X  coordinate  than  when  the  aircraft  has  a  positive  X  coordinate.  An 
approximation  to  the  drogue  yo-yo  is  the  difference  between  the  longest 
and  the  shortest  straight  line  distances  from  the  aircraft  to  the  junc¬ 
tion  measured  in  the  X-Z  projection.  The  longest  distance  occurs  when 
the  aircraft  :'.o  located  at  the  position  on  its  orbit  coincident  with 
the  negative  X  axis,  and  the  shortest  distance  occurs  when  the  aircraft 
is  positioned  at  the  positive  X  axis.  Applying  this  approx.  ..ate  method 
of  determining  the  drogue  yo-yo,  it  is  difficult  to  locate  the  junction 
separating  the  uppermost  and  lowermost  sections  of  the  cable  without 
first  viewing  a  time  exposure  of  the  dynamic  motion  of  the  cable.  Of 
course,  the  junction  must  be  within  the  envelope  formed  by  the  sus¬ 
pended  cables.  It  is  sufficient  to  assume  that  the  junction  is  located 
on  a  curved  line  parallel  to,  and  midway  between  the  suspended  cables. 

Now  the  determination  of  the  yo-yo  is  reduced  to  positioning  the  junction 
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on  the  curved  line.  Figure  17  is  a  graphical  representation  of  the 
approximate  method.  3ased  on  the  experience  gained  from  the  time  ex¬ 
posure  observations  a  good  approximation  of  the  junction  point  is  some¬ 
where  between  a  quarter  and  a  half  way  down  the  suspended  cable  shape. 


I  The  cable  shape  inferred  here  is  the  dashed  line  which  is  the  computer 
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predicted  suspended  cable  shape  based  on  the  parameters  listed  on  page 
52.  On  Figure  17,  point  C  was  used  for  the  yo-yo  calculation.  To  indi¬ 
cate  the  freedom  involved  in  positioning,  the  junction  point  on  the  sus¬ 
pended  cable  shape,  points  C'  and  C"  were  also  used  which  resulted  in 
less  than  3%  variation.  In  order  to  calculate  the  drogue's  average  des¬ 
cending  and  ascending  speeds,  one  divides  the  yo-yo  by  the  length  of 
time  the  aircraft  takes  to  fly  the  upwind  and  downwind  portions  of  the 
orbit.  By  programming  the  aircraft  to  orbit  with  the  velocities  found 
in  Section  2.2,  these  lengths  of  time  are  easily  determined.  In  apply¬ 
ing  this  approximate  method  to  ascertain  the  yo-yo  and  average  speeds 
of  the  drogue,  ell  necessary  calculations  can  be  made  on  a  programmable 
desk  calculator,  since  finding  the  shape  of  the  suspended  cables  in¬ 
volves  static  equations  of  equilibrium  ai.d  is  trivial  compared  to  the 
integration  of  the  dynamic  equations. 


Returning  to  Figure  20  and  using  the  approximate  method  described 
above,  the  drogue  yo-yo  is  calculated  as  1,030  ft  when  the  critical 
junction  is  located  1,600  ft,  or  one-fifth  of  the  entire  cable  length 
from  the  orbit.  Thic  distance  is  measured  along  the  curved  line  paral¬ 
lel  to,  and  midway  between  the  two  suspended  cables.  The  1,600  ft 
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distance  is  selected  since  it  gives  accurate  values  of  yo-yo  for  each 
of  CASES  3-6.  Using  this  value,  the  junction  appears  to  be  within  the 
uppermost  section  of  the  cable.  The  yo-yo  calculated  by  integrating 
the  equations  of  motion  for  the  system  is  1,247  ft.  The  approximate 
value  of  yo-yo  is  13.6%  less  than  the  exact  value. 

CASE  4 

Applying  the  approximate  method  to  Figure  21,  the  yo-yo  turns  out 
to  be  680  ft  when  there  are  1,600  ft  separating  the  critical  junction 
from  the  aircraft  orbit.  The  exact  yo-yo  is  682  ft,  so  the  approximate 
yo-yo  is  .3%  less  than  the  exact  value. 

Note  in  Figure  21  that  there  is  no  flaring  at  the  bottom  of  the 
funnel-shaped  surface  in  a  15  KTAS  wind. 

CASE  5 

The  same  approximate  method  is  used  to  find  the  yo-yo  when  flying 
an  elliptical  orbit.  Since  the  X-Z  projections  of  both  the  orbit  and 
the  envelope  seen  in  Figure  22  are  the  same  as  in  Figure  21  of  CASE  4, 
then  680  ft  of  yo-yo  will  again  be  calculated  when  1,600  ft  of  cable 
seo.irate  the  aircraft  orbit  frcm  the  critical  junction.  This  is  7.6% 
less  than  the  exact  value  of  736  ft. 

CASE  6 

For  the  elliptical  orbit  in  Figure  23,  the  approximate  method  pre¬ 
dicts  1,000  ft  of  yo-yo  when  1,600  ft  of  cable  separate  the  aircraft 
from  the  critical  junction.  This  is  1.8%  greater  than  the  exact  value 


of  yo-yo  which  is  990  ft. 


CASE  7 

In  Figure  24  there  is  noticeable  flaring  at  the  bottom  of  the 
surface  swept  out  by  the  4,000  ft  cable,  but  no  flaring  at  the  bottom 
of  the  surface  generated  by  the  8,000  ft  cable.  Further  analysis  of 
these  surfaces  in  connection  with  the  approximate  method  of  determining 
yo-yo  explains  why  the  yo-yo  is  nearly  the  same  for  both  systems.  The 
vertical  separation  between  the  aircraft  orbit  and  the  critical  junction 
is  greater  for  the  8,000  ft  cable.  If  the  junction  had  the  same  X 

coordinate  for  both  cables,  the  yo-yo  would  be  less  for  the  8,000  ft 

cable.  However,  the  4,000  ft  cable  is  not  blown  as  far  downwind,  and 
its  yo-yo  is  slightly  less. 

CASE  8 

Increasing  the  weight  of  the  drogue  from  25  to  100  lbs  causes  the 
flaring  at  the  bottom  of  the  surface  to  become  appreciable,  as  seen  in 
Figure  25.  The  approximate  method  of  determining  yo-yo  explains  why 
the  heavier  system  has  the  least  amount  of  yo-yo.  The  vertical  separa¬ 
tion  between  the  aircraft  orbit  and  the  critical  junction  is  greater 
for  the  system  employing  a  100  lb  drogue.  This  would  cause  the  yo-yo 

to  be  less  for  the  heavier  system,  if  the  junctions  had  the  same  X  co¬ 

ordinate.  Since  the  heavier  system  does  not  blow  as  far  downwind,  the 
yo-yo  is  even  further  reduced. 

CASE  9 

In  Figure  26  the  surface  generated  by  the  8,000  ft  r*able  towed  in 
a  nonconstant  altitude  orbit  is  most  difficult  to  distinguish  from  the 
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h  »  18,400  ft 

Cable  density  ■  234.7  lb/ft^ 

Cable  diameter  ■  .17  in 

6  2 

Cable  elastic  modulus  ■  10x10  lb/in 
Cable  length  *  16,000  ft 
Drogue  weight  »  250  lb 
Drogue  radius  (spherical)  =»  2.5  ft 

The  yo-yo  calculated  by  integrating  the  equations  of  motion  with  N=20 
equals  3,825  ft.  Figure  17  graphically  demonstrates  the  utilization 
of  the  approximate  method.  The  resulting  yo-yo  is  4,000  ft,  or  4.4% 
greater  than  the  exact  value  of  yo-yo.  This  system  has  pointed  out  the 
versatility  of  the  approximate  method  of  calculating  the  drogue  yo-yo. 
It  can  be  applied  with  adequate  accuracy  to  almost  any  orbital  trailing 
wire  system. 

The  surfaces  generated  by  the  towed  cables  flown  in  prevailing 
winds,  as  seen  in  Figures  20-26,  illustrate  the  inaccuracy  of  a  simple 
model  of  the  orbiting  trailing  wire  system  that  was  proposed  in  the 
University  of  Dayton  Research  Institute  Technical  Report  UDRI-TR-70-31 
of  July  1970.  The  model  is  described  in  the  Technical  Report  as 
follows : 

"The  model  consists  of  assuming  that  centers  of  the 
circles  made  by  each  segment  of  the  cable  assume  a 
shape  similar  to  that  of  a  motionless  cable  suspended 
from  the  center  of  the  aircraft  orbit.  The  surface 
in  which  the  cable  moves  is  then  constructed  around 
this  line." 

(Crouch,  Bauer,  Kahle,  1970,  p.  28) 
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CONCLUSIONS 


a.  A  limped  mass  model  of  trailing  wires  can  adequately  predict 
the  dynamics  with  a  relatively  few  mass  points  and  with  a  reasonable 
expenditure  of  computer  time. 

b.  Both  cross  winds  and  orbital  perturbations  can  cause  the  yo-yo 
phenomenon  in  cables. 

c.  In  general,  long  cables  towed  from  nonconstant  altitude  cir¬ 
cular  orbits  from  slowly  flying,  steeply  banked  aircraft  with  moderate 
weight  drogues  will  give  the  most  desirable  cable  characteristics. 

d.  Static  analysis  of  Figure  17  adequately  predicts  the  amount 
of  yo-yo  that  will  occur. 
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Figure  8a 
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Figure  9a 


Figure  9b 
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Figure  1  la 


Figure  1  lb 
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Figure  13a 


Figure  13b 
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TILTED  ELLIPTICAL  FLIGHT  PATH 
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APPENDIX  I 


The  cable  is  divided  into  n  mass  points  connected  by  straight, 
linear  elastic,  massless  segments.  If  an  adequate  number  of  mass  points 
are  taken,  the  model  accurately  predicts  the  response  of  the  actual 
cable.  Each  mass  point  has  three  degrees  of  freedom  and  is  located  in 
the  x,  y,  z  coordinate  system  shown  in  Figure  27. 

The  equations  of  motion  can  be  derived  using  Lagrange’s  Equation 
and  treating  the  aerodynamic  forces  as  generalized  forces.  The  kinetic 
energy  of  the  cable  system  is  given  by 
n 


m. 


V1  i  ,*2  .  *2  ,  * 2. 

2-1  T  (xi  +  yi  +  zi>  • 


i«l 


the  potential  energy  is  given  by 


(1) 


v  -  E  r  (ei "  so )2  +  E  v, 

i=l  1  i=l 

where  m^  *  mass  of  iC^  point 

k^  =  spring  constant  for  segment 


(2) 


e^  =  stretched  spring  length 


sq  =  unstretched  spring  length 


g  =  32.3  ft/sec  . 

For  any  given  mass  point,  i,  the  two  connecting  segments  e^  and 
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are  given  by 


In  addition 

3e<  (^x-rO 


i  _  "‘i+l  Hi'  ,  3ei-l  _  (qi‘qi-l) 


**4 


and 


3*4 


'i-1 


where  again  the  q^’s  are  x^,  y^  or  z^. 


The  last  step  is  to  calculate  the  Q^  forces  due  to  lift  and 
drag  on  each  cable  segment.  The  wind  velocity  vector,  V^,  is  due  to 


the  cable  motion  and  can  be  expressed  as 

^  ^  i  -  y±  j  -  ^  ic  . 

the  vectors  along  the  e^  and  e^_x  segments  can  be  expressed  as 

*i =  (xi+i  -  V  1  +  (yi+rV  5 +  (zi+r2i)  k 

and 

A  A  A 

ei-l  =  ^xi~Xi-l^  1  +  ^yi_yi-l^  $  +  ^Zi-Zi-1^  k  * 

From  the  dot  product,  the  angles  of  attach,  g^  and  g  ^  can  be 

found  to  be  given  by 


(8) 


(9) 


S  =  cos 


-17  Vi'ei 


Mini)  (10)' and  3i-i '  c°rl  ngft  •  (11) 


From  the  cross  products  one  can  find  the  vectors  in  the  direction 
of  lift  (i.e.  perpendicular  to  and  in  the  plane  formed  by  and  e^ 

or  e^_n).  Letting  and  ri^_^  be  unit  vectors  in  the  lift  direction  we 
get 


(V.xe)xV  .  (V.xe  )x  V 

n.  -  '  Z~Z - Z  (12)  and  =>  -3-“ - ±- 

1  |(Vixe1)xVij  |VIxei_1)xVij 


(13) 
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